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Abstract 

An external magnetic field has been applied in laterally coupled dots (QDs) and we have studied 
the QD properties related to charge decoherence. The significance of the applied magnetic field to 
the suppression of electron-phonon relaxation and dephasing rates has been explored. The coupled 
QDs have been studied by varing the magnetic field and the interdot distance as other system 
parameters. Our numerical results show that the electron scattering rates are strongly dependent 
on the applied external magnetic field and the details of the double QD configuration. 
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I. INTRODUCTION 



Relaxation and decoherence properties in electronic devices such as low-dimentional struc- 
tures and specifically in quantum dots (QDs) have attracted the interest of experimental and 
theoretical studies. Fabricated coupled QDs have been suggested as candidates for quantm 
bits (qubits) where the states of a single trapped electron within the coupled QDs can per- 
form the two states of the qubit.-'^'^'^i^'^ The investigated devices are mainly characterized 
by two important decoherence channels. The first one is the coulomb interaction to the 
background charge fluctuation which is an extrinsic decoherence and the second channel, 
an intrinsic decoherence, due to electron-phonon interaction. The theoretical studies of the 
above mentioned intrinsic decoherence due to single elctron-phonon coupling are mainly 
related to the single phonon emission process. Although the electron relaxation through 
multi phonon processes is an important parameter at high temperature, at qubits operation 
temprature (~ IT) these processes have not strong effects.- 

It is worth mentioning that the last decade a vast of research has been published on 
double QDs dephasing^. In our previous work^, we have reported the results of charge 
decoherence due to relaxation and dephasing rates caused by electron-phonon interaction. 
In the present investigation, we report the significance of longitudinal acoustic (LA) phonons 
(including the deformation and piezoelectric interactions) and LO phonons on relaxation 
and dephasing rates in laterally coupled QDs under the existence of magnetic fields. The 
existence of an external magnetic field is of special importance for the electron states^ due 
to the decrease of electron energy splitting as the magnetic field increases. The scattering 
rates have been found to depend strongly on electron confinement and the interdot distance. 
Acoustic phonons are only used in the calculations of the relaxation rate due to the small 
electron energy splitting (<lmeV). For the dephasing rates, both the acoustic and optical 
phonons are considered. 

The paper is organized as follows. In section II, we describe the electron wavefunctions 
under the present of an external magnetic field and the models for the acoustical and optical 
phonon modes. In section III, we outline the scattering theory of an electron which is 
scattered to an energetically smaller state by the emission of acoustical phonons. The 
dephasing rates have been also studied by considering the emission of acoustical and optical 
phonons. Section IV is devoted to the numerical results and interpretation of suppression 
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of decoherence due to magnetic field. Lastly, section V presents a summary of the present 
work. 



II. THE MODELS OF CHARGE QUBIT AND PHONONS 
A. Electron wavefunctions 

We consider a single electron within two coupled identical QDs with an interdot distance 
la. The lateral electron motion is decoupled from the one along the QW growth.^- The 
Hamiltonian which describes the single electron motion which is confined in laterally coupled 
QDs is given by 

n = n\\^n. (i) 

where the QW growth is denoted by the subscript "2;" and the lateral direction by "||". 

The external magnetic field enters the Hamiltonian via a magnetic vector potential A. 
By choosing the symmetric gauge K='B (— + xey) /2 then the magnetic field points to 
z direction and it is given by B = V x A = iSe^. The lateral confinement is assumed to be 
parabolic for a single QD, therefore, the Hamiltonian operator for the lateral directions has 
been considered as^ 

Till = ^ + -m*wV[? - -cJcL^ (2) 
" 2m* 2 II 2 ^ ' 



where the operator of the z component of the angular momentum is given by 



d d' 
ox oy 



(3) 



with p the quantum mechanical operator of momentum, ujq is a parameter (in this case 
frequency) describing the strength of the confinement in x-y plane, Uc = Be/m* and uj"^ = 
u!q + {uJc/2)'^ . The electron wavefunction can be separated to the following envelope functions, 

=^ll (i'll)^-(^) (4) 

In the case of a single QD, the one electron wavefunction can be given in terms of the 
principal quantum number n (n = 0, 1, 2, ...) and the angular momentum quantum number 
m (m = 0, ±1, ±2, ...) as 



(n,m) ,~ n\ _ ^- ~\m\ -^^2 imO r^\m\ ( ~2\ (r\ 



where £[^1 (p^) are the Laguerre polynomials, and p is a scaled radius, p — with 



/ = yh/m*uJo. The eigenvalues of the single particle are given by 

Km = (2n + |m| + 1) huQ (6) 

By using the Heaviside step function 6, the Hamiltonian, along the QW growth, takes the 
form 

= + VoQ {\z\ - L,) (7) 

2 m*[z) 

where m*{z) is the electron effective mass and Vq the offset between the band edges well and 
barrier. The wavefunction along the above mentioned direction has been considered as the 
wavefunction of an infinite QW (Vq oo). Here, it has been used only the ground state 
wavefunction along the QW growth due to the strong confinement along this direction. The 
ground state wavefunction is given by ipz (z) — Acos{t:z/2Lz) where ^ is a coefficient to 
be determined by normalization and 2L^ is the size of the QW. 

The wavefunction of a single electron which is confined in a 2D QD and is coupled in one 
dimension of the x-y plane, can be formed by superposition of two uncoupled QDs which 
are sited "left" and "right" of the the origin of the frame of reference and are separated by 
an inter-dot distance 2a. The external confining potential which it is used is given by 

Vc—-m*UQrriin{{x — af-\-'lf'., (x + a)^ + y^} (8) 

The single electron wavefunction for the parallel plane can be given by: 

*ii) = E<^^K^)+^^K«) (9) 

k 

and the total wavefunction of the system of coupled QDs as described above is 

*(r) = *||(r||)v^,(^) (10) 

The wavefunctions in the parallel plane of the coupled QDs system are calculated numerically 
by direct diagonalization. 



B. Acoustic and optical phonons 



The electrons, in polar semiconductors couple to acoustical and optical phonons. The 
small electron energy splitting does not permit any electron transition in the charge qubit 
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via the emission of optical phonons. On the other hand, the acoustical phonons contribute to 
the relaxation rates due to the small phonon energies. Here, we calculate the decohence rates 
which are caused due deformation potential and piezoelectric acoustic phonon interaction 
by considering only longitudinal phonons. It follows the Hamiltonian which describes these 
interactions^^: 

^=i:(^)"'A<wM)K+"i,). (11) 

where cUq is the frequency of the phonon mode with wavevector q, pm is the mass density of 
the host material, V is the volume of the sample, Oq and phonon annihilation and 

creation operators, and p(q) is the electron density operator. The term A^(q) is given by 

M{ci) = D\ci\ + iMx{q) (12) 

The first term of the above equation represents the deformation potential interaction with 
deformation potential D and the second part, which is imaginary, is the piezoelectric 
interaction.— For zincblende crystals (e.g. GaAs), the term M.x{<i) can get the form^^ 

Alf (q) = 2e ei4 {q^Qyiz + Wzix + qxQziy) (13) 

where e is the electronic charge, ei4 is the piezoelectric constant, and ^ is the unit polarization 
vector. 

The fact that the energy difference between the electron states in coupled quantum dots 
is quite small (a few meV) does not allow optical phonon transitions due to the conservation 
of energy (optical phonon energy is ~ 36 meV in GaAs). The optical phonons play a role in 
dephasing rates as we demonstrate in the next section. By using the bulk phonon approxi- 
mation and neglecting the interface phonon modes^^ii^, the electron-phonon interaction due 
to LO phonons is thus given by^ 

Hop = E ^^(q) («q + «-q) (1^) 
q qW 



= 27ce^hujLo ( — - - ) (15) 



and 

= 27re^hujLo ( 

V Coo 

where ulo is the longitudinal optical frequency, and e^o are the static and high frequency 
dielectric constant of the host material. 

Having described the electronic states in coupled QDs and the relevant types of phonon, 
we are now ready to calculate the relaxation and dephasing rates. The next section is 
devoted to Fermi's golden rule and dephasing rates. 



III. THEORY OF RELAXATION AND DEPHASING RATES 



The relaxation rate between an initial 



'^^^^) and a final state "^^^n associated with 



phonon emission (+) (or absorption (-)) is determined by Fermi's golden rule: 

T 



r = ^^|(^(^)(r)|/7'"*k(^)(r))|'5(E^-E,±E, 



q 



NB{E^,Tiat) + l±l) (16) 

where the labels T and 'F' denote the initial and final electron states respectively. A^^ is the 
Bose- Einstein distribution function for phonons with lattice temperature Tiat- It is worth 
mentioning that in our calculations we assumed Tiat = and the phonon absorption can be 
neglected. 

It is also worth mentioning that relaxation is not the only way charge qubits can be deco- 
hered. If the energy difference between the two charge states fluctuates, phase information 
will get lost and decoherence occurs. The density operator of an electron in a boson bath is 
given ini^ii^ 

^ Poo(O) poi(0)e-^^(^*)+-^*/^ 

^Pio(0)e-^'(^*)-^*/^ pn(0) 
where e is the energy splitting between the electron energy levels. In short, pure dephasing 
cause a decay in the off-diagonal element of the density matrix for the two-level system that 
makes up the charge qubitr^'^^: 

Poi(t)~Poi(0)e-^'W, (18) 
where the exponent function B'^{t) is deflned by 

For acoustic phonons, we choose frequencies cUq = qCg for the relevant branches, while for 
longitudinal optical phonons, we choose c<jq = ujlo- The coupling constants g{q) due to 
deformation potential, piezoelectric and optical phonons are respectively given by 



Pit) 



(17) 



^piezo(q) = -Mr(q) /^^(q) ' (21) 

^poiar(q) = ^^X(q) , (22) 
qVV 
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where X(q) is given by 



(23) 



the symbols (±) refer to the two states for the double dot charge qubit. The matrix integrals 
in this study are carried out using the Monte-Carlo algorithms. 



IV. RESULTS AND DISCUSSIONS 



We firstly calculate the relaxation rates as a function of an external magnetic field B, for 
an electron which scatters from the first excited to ground state associated by the emission 
of acoustical phonons. In all our calculations, the quantum well width is fixed to the value 
of 2Lz = 10 nm and the confinement lengths in the x an y directions are 0.5 fim. 

As it appears in Fig. [H for small magnetic field the relaxation rates due to deformation 
interaction is larger than the one due to piezoelectric interaction. As the field B increases 
the piezoelectric coupling becomes the dominate contributor due to the different wavevec- 
tor dependence in the deformation and piezoelectric matrix elements. For the deformation 
potential, the dependence of matrix elements on wavevector is related to y/q while for the 
piezoelectric coupling is related to 1/^^. The relaxation rates for small magnetic field in- 
crease up to a maximum value and afterwards decrease. This resonance reflects the existence 
of large electron wavefunctions of the charge qubit (for B ~3 Tesla). 

Fig. [2] presents the electron relaxation rates as a function the half the interdot distance for 
a fixed B = 3 Tesla and hu = 3 meV. As the interdot distance increases the rates decrease 
due to the small energy splitting which results small phonon density of states. Thus the 
relaxation rates become small as the interdot distance increases. For a close to 21.5 nm, 
it is obvious that the rates get a maximum value due to large wavefunctions. The different 
dependence of matrix elements on the phonon wavefuction for the deformation potential and 
piezoelectric interaction results the different behavior of the relaxation rates for the above 
mentioned interactions. 

The dependence of the relaxation rates on the electron confinement strength is shown 
in Fig. [3l The rates increase as the energy splitting between the first excited and ground 
state increases. When the electron strength of the electron confinement reaches the value 
huj = 6.5 meV, the relaxation rates get a maximum value. Increasing the electron confine- 
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ment strength, the energy sphtting becomes small and the relaxation rates decrease as a 
consequence of the energy splitting degradation. The dependence of relaxation rates on the 
energy splitting is shown in inset [21 Which reflects the fact that for a given energy splitting, 
correspond two values of the strength of the confinement.^ As a result there are two values 
of relaxation rates for each energy splitting. The two different types of electron-phonon in- 
teractions produce different contribution to the total relaxation rates and can be interpreted 
in the same manner as in Fig. [TJ 

The second part of our investigation of decoherence in charge qubits is the evaluation of 
dephasing factor and its dependence on an external magnetic field. Here, we calculate the 
dephasing effects from both acoustic (deformation, piezoelectric) and optical phonons. Fig. 
m shows the temporal dependence for two different values of an external magnetic field. The 
curves in Fig. H] rapidly increase for the first 10 ps and and for later times they saturate. 
As a result -B^(t) depends only very slowly on time after 100 ps. The interaction between 
the qubit electron and the acoustic phonon bath causes the fast increasing of dephasing in a 
period of time less than 100 ps. Mathematically, the very fast time dependence of dephasing 
is due to the trigonometric dependence on phonon frequencies and time [see Ref. '9(]. For 
larger time is practically flat and can be considered constant after 100 ps. A constant 
dephasing factor will not produce a decaying signal in terms of, for example, oscillations in 
electrons. Instead, it simply reduces the contrast in the charge oscillation. This can be seen 
easily from Eq. [171 The presence of a constant exp{—B^) ~ exp(— 0.05) simply reduces the 
magnitude of poi by a constant factor of 0.05, which is not a particularly large suppression 
(though significant in terms of fault tolerant quantum computing). 

The temporal behavior of the dephasing factor has an interesting feature for zero and a 
none zero magnetic field. B'^{t) decreases as B increases (Fig. W^,h) due to the dependence of 
the matrix elements on the initial and final wavefunctions under the presence of an external 
magnetic field. Fig. [5] shows the dephasing rates as a function of an external magnetic 
field. As the field B increases the quantity X(q) in Eq. (1231) decreases as a consequence of the 
smaller differences in the matrix elements involved in Eq. (l23i) . For large magnetic fields 
the rates go to zero. 

Finally, we calculate the dephasing rates as a function of the half the interdot distance 
for a fixed external magnetic field and electron confinement strength. As in the case of the 
relaxation rates, the dephasing rates decrease by increasing the interdot distance due to the 
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fact that when the two QDS are well separated then the integral difference in Eg. (1231) goes 
down quickly. Consequently, the dephasing factor B^{t) undergoes suppression by increasing 
the separation distance. 

V. CONCLUSIONS 

We have studied the phonon-induced single electron relaxation and dephasing rates in 
laterally coupled QDs with the presence of an external magnetic field. The relaxation 
and dephasing rates have been calculated for different system parameters such as interdot 
distance, strength of electron confinement and magnetic field. In the case of zero magnetic 
field^., the rates could be enhanced for some double dot configurations. This enhancement 
of the rates could be easily suppressed by the existence of an external magnetic field. Our 
results show that the magnetic field is of crucial importance in the study the decoherence 
in charge qubits due to the suppression of electron relaxation and dephasing rates. 
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FIG. 1: Relaxation rates of an electron, which scatters from the first excited state to ground state, 
versus the external magnetic field. The confinement strength is Tiuiq = 3 meV, half the interdot 
distance a = 20 nm and QW width 2Lz = 10 nm.. 



4x10 




a (nm) 

FIG. 2: Relaxation rates of an electron versus half the interdot distance a. The confinement 
strength is hujo = 3 meV, the magnetic field B = 3 Tesla and QW width = 10 nm. 
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FIG. 3: Relaxation rates of an electron versus the strength of the confinement Huq. Electron 
relaxation rates as a function of the energy splitting between the first excited state and the ground 
state is given in the inset. The half the interdot distance is a = 20 nm, magnetic field B = 3 Tesla 
and QW width 2L^ = 10 nm. 
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FIG. 4: Dephasing factor B'^{t) as a function of time t. The dephasing factor due to deformation 
potential, piezoelectric interaction, LO and the total rates are represented by dashed, dotted, 
dashed-dotted and straight line respectively. The strength of the confinement is Tiujq = 3 meV, half 
the interdot distance is a = 20 nm, with a) H = Tesla and b) H = 3 Tesla. 
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FIG. 5: Dephasing factor B'^{t) as a function of the magnetic field. The time t is chosen to 
be 60 ps. The sohd hne represents dephasing rates due electron-acoustic phonon interactions. 
The dephasing rates due to deformation potential, piezoelectric interaction and polar interaction 

with optical phonons are represented by dashed, dotted, and dash-dotted line respectively. The 
strength of the confinement is htoo = 3 meV, half the interdot distance a = 20 nm and QW width 
2Lz = 10 nm. 
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FIG. 6: Dephasing factor -B^(i) as a function of the half interdot distance a. The time t is chosen to 
be 60 ps. The solid line represents dephasing rates due electron-acoustic phonon inter actions. The 
dephasing rates due to deformation potential, piezoelectric interaction, and polar interaction with 
optical phonons are represented by dashed, dotted, and dash-dotted line respectively. The strength 
of the confinement is hiOo = 3 meV, B = 3 Tesla and QW width 2Lz = 10 nm. 
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